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Abstract 

2^ ■ Using the data for 24.5 x 10^ il^C^S) produced in e~^e~ annihilations at ^/s = 3686 MeV at 

• ! the CESR-c e'^e~ colhder and 8.6 x 10^ J/ip produced in the decay ip{2S) — t- vr+vr^J/^, the 

^^ I branching fractions for 'ip{2S) and J/ip decays to pairs of pseudoscalar mesons, vr+vr^, K^K^ , and 

KsKi, have been measured using the CLEO-c detector. We obtain branching fractions B{ip{2S) — )• 
vr+vr-) = (7.6 ± 2.5 ± 0.6) x 10-^ 13{'^j{2S) -^ K+R-) = (74.8 ± 2.3 ± 3.9) x 10-^ B{i){2S) -^ 
KsKl) = (52.8 ± 2.5 ± 3.4) x 10-^ and i3( J/V' ^ vr+vr-) = (1.47 ± 0.13 ± 0.13) x lO^^, B{J/tp -^ 
K+K-) = (2.86 ± 0.09 ± 0.19) x 10-^ S(J/V' ^ i^^i^i) = (2.62 ± 0.15 ± 0.14) x 10-^ where 
the first errors are statistical and the second errors are systematic. The phase differences between 
the amplitudes for electromagnetic and strong decays of ip{2S) and J /ip to '' pseudoscalar pairs 
are determined by a Monte Carlo method to be 5{'il:{2S)pp = (110. 5+9 g )° and 5{J /'ip)pp = 
(73.5t|g)°. The difference between the two is AJ = 5{ip{2S))pp - 6{J/ip)pp = {37.0l:\o.l)°- 



PACS numbers: 13.25.Gv,13.66.Bc,14.40.Be,14.40.Pq 



I. INTRODUCTION 

Interest in final state interaction (FSI) phases originally arose from CP violation in K 
decays and B decays. However, recently interest has focused on the suggestion that large FSI 
phases are a general feature of strong decays. The electromagnetic and strong decays of the 
vector states of charmonium, J /il) and ■?/'(2S'), offer good testing ground for this possibility. In 
a series of papers Suzuki [1-3] has studied FSI phase differences between the electromagnetic 
amplitude A^ and the three-gluon strong amplitude Aggg. In the decays of J/ip to 1~0~ 
vector-pseudoscalar (VP) pairs Suzuki obtained 6{J/iIj)vp = 80° [1], and for the 0~0~ 
pseudoscalar-pseudoscalar (PP) pairs he obtained S{J/ip)pp = (89.6 ± 9.9)° [2]. Rosner [4] 
confirmed Suzuki's results, obtaining 5{J/ip)vp = (76^^q)°, and 5{J/4')pp = (89 ± 10)°. 
Further, Gerard and Weyers [5] have argued that these differences are manifestations of 
what they call "universal incoherence", i.e., 90° phase difference between electromagnetic 
and every exclusive annihilation decay channel of J/ip and iIj{2S). In order to arrive at 
a deeper understanding of the origin of large phase differences it is therefore necessary to 
examine if what has been observed for J/ip decays persists in the corresponding decays of 
iIj{2S). As Suzuki has noted [3], this is particularly important in the context of the curious 
suppression of the ratio r[ilj{2S) — ?■ VP]/r[J/ilj — )> VP] (particularly notable for pn decays). 

One of the best places to study the phase difference between electromagnetic and strong 
decays oiilj{2S) and J/ip is in their decays to pseudoscalar pairs, tt+tt^, K^K~, and KsK^. 
This is because the three decays sample the interactions in quite different ways. The tt+tt" 
decay is essentially purely electromagnetic, with strong decay being forbidden by isospin 
invariance, the KgKi decay is essentially purely strong and due only to SU(3) violation, 
and the K^K~ decay can proceed through both the electromagnetic and strong interactions. 

A particularly simple and transparent determination of the relative phase angle 5{il)) can 
be made by measuring it as the interior angle of the triangle in the complex plane with the 
amplitudes for the decays to tt^tt", K^K~ , and KsK^ as its three sides. This is illustrated 
in Fig. 1 for the J /ip and ip{2S) decays. 

With the neglect of the small effect of SU(3)-breaking and interference between resonance 
and continuum amplitudes, the relative phase angle 6{iIj) is given by 

V \2V^iKsKL)xpxB{n+7r-)\ i ' 




A„„„ Qc [Br(J/v-*KsKJ]i'2 




A^, oc [Br(>|;(2S)^KsKL)] 



FIG. 1. Triangle representations of the amplitudes in J /ip — t- PP and ip{2S) — )■ PP decays. The 
relative phase angle between the strong and electromagnetic amplitudes is 6{ip). 



where the phase space correction factor p = {pk/PttY'i p(V^(25')) = 0.902, and p{J/ip) = 
0.862. Thus it is required to measure the three branching fractions, B{ip{2S), J/ip — )■ 
TT^TT^, K^K~ , and KsKl), which are proportional to the squares of the respective ampli- 
tudes. 

Previous measurements of 6{iIj{2S))pp were made by BES [6, 7] and CLEO [8]. Their 
results, recalculated to correspond to the internal phase angle defined by Eq. (1), are 
6{i/j{2S))pp = (91 ± 35)° (BES), and 6{ip{2S))pp = (87 ± 20)° (CLEO). The results for 
S{J/tp)pp mentioned earlier were obtained by Suzuki and Rosner using the PDG 1998 [9] 
evaluation of the branching fractions measured in 1985 by Mark III [10] in the decay of their 
sample of 2.7 x 10^ J/ip produced directly in e~^e~ annihilations at ^/s = M{J /ip). In this 
paper we present much more accurate results for branching fractions and 5{il){2S))pp using 
the CLEO-c data of 24.5 x 10^ il){2S)^ eight times larger than that used in the previous 
CLEO measurement [8], and for branching fractions and S{J/iIj)pp using the data set of 
8.6 X 10^ J/ip tagged by tt+tt" recoils in the decay ip{2S) — )■ -k^-k^ J/ip. In all cases large 
improvement in precision over previous measurements is obtained. 

The data used in this analysis were collected at the CESR e+e^ storage ring using the 
CLEO-c detector [11]. The detector has a cylindrically symmetric configuration, and it 
provides 93% coverage of solid angle for charged and neutral particle identification. The 
detector components important for the present analysis are the vertex drift chamber, the 
main drift chamber (DR), the Ring Imaging Cherenkov (RICH) counter, and the Csl crystal 
calorimeter (CC). 



II. EVENT SELECTION 

Event selection, particle identification, and branching fraction determination in the 
present paper follow closely those used by us in our published paper on the determination 

oi6{^{2S)) [8]. 

The events for the three decay modes, iIj{2S), J/ip — )> tt+tt^, K^K' , Ks{^ 7i^7i^)KL{not detected) 
are required to have two charged particles for ip{2S) decays, and four charged particles for 
J/ip decays, and zero net charge. The event selection criteria are identical to those in our 
earlier paper [8]. To recapitulate, all charged particles are required to meet the standard 
criteria for track quality. The tt+tt^, K^K^, and tt+tt^ from Ks decays (with vertex dis- 
placed by > 5 mm) are accepted in the regions | cos^| < 0.75, 0.80 and 0.93, respectively. 
The invariant mass of the tt+tt^ from Ks decay is required to be within ±10 MeV of 
M{Ks) = 497.61 MeV. 

Particle identification is done by combining dE/dx information from the drift cham- 
ber (DR) and the likelihood information from the RICH detector for the particle species 
2, j = p, K, 71, fi, and e. The variable for the dE/dx information is Si = [{dE/dx)j^casurcd — 
{dE / dx)cxpcctcd]/ o'ldE / dx)raca.suvcd, and for the RICH information it is the likelihood func- 
tion, -21ogL. To distinguish between two particle species a joint x^ function, Ax^(i, j) = 
— 2(logLj — logLj) + {Sf — S'j) is constructed. Charged kaons in K'^K~ decays are distin- 
guished from protons, pions, and leptons by requiring I\')^{K,'pJT\ j ^ije) < —9. Looser 
criteria are used for pions in tt+tt^ decay and tt^tt^ daughters from Ks in KsKl de- 
cay, Ax^iT^, e/K/p) < 0. In addition, electrons are rejected in all decays by requiring 
E{CC)/p < 0.9. All these requirements are identical to those in Ref. [8]. 

In Ref. [12] a detailed study was made to distinguish pions from the much more prolific 
yield of muons. It was determined that the energy deposited in the central calorimeter by 
pions due to their hadronic interactions provides a very efficient means of distinguishing them 
from muons which deposit much smaller energy due only to dE/dx. For ip{2S) — )> tt+tt" 
decay, it was determined that requiring every pion to deposit E{CC) > 0.42 GeV reduced 
muon contamination to ^ 1% level. This requirement was used in Ref. [8], and we impose 
it also in the present analysis for the channel ip{2S) — > tt+tt^. For J/ip — )> tt+tt", the pion 
yield is much larger, and the corresponding requirement is determined to be E{CC) > 0.35 
GeV. For KsKl decay an explicit n^ veto was made, as in Ref. [8]. 



III. DETERMINATION OF B(ipi2S) -^ tv+tt', K+K', KsKl), AND d(^(25)) 

For ip{2S) — 7> TT+TT^ and K^K~, it was required that the total momentum |Ep| be less 
than 75 MeV. As expected, this requirement removes most of the J/ip, Xcj peaks and the 
background which is present without it in the event distributions plotted as a function of 
X{h) = {E{h^) + E{h'))/^/s. The resulting X{h) distributions are shown in Figs. 2(a,b) 
in the extended region of X{h). For ip{2S) — )■ KsKl with only Kg observed no such 
total momentum cut can be imposed, and the different criteria developed in Ref. [8] were 
implemented to take account of the unobserved Kl. 

For ip{2S) — )■ KsKl, the direction of K^ is inferred from that of the observed Ks- We 
require that there be no shower associated with neutrals closest to this direction with energy 
> 1.5 GeV. Further, we require that in a cone of 0.35 radians around the Kl direction there 
be no more than one shower with i^in > 100 MeV. We require that outside this cone there be 
no single shower with i^out > 100 MeV and the sum of all showers Si^out < 300 MeV. These 
selections remove events with neutral particles other than K^ accompanying the detected 
Ks- As shown in Fig. 2 (c), the remaining background is featureless, and very small in the 
signal region, X{Ks) ~ 1.0. 

In Figs. 2 (d,e,f) we show the X{h) distributions of Figs. 2 (a,b,c) in detail in the smaller 
region of X{h) in which we fit the data with MC-generated peak shapes and linear back- 
grounds. The fit results are presented in Table I. The observed peak counts, A^(obs), are 
70.8 ± 8.8, 1431.3 ± 39.4, and 478.0 ± 23.0 for tt+tt-, K+K", and KsKl, respectively. 

The MC signal events were generated assuming sin^ 9 angular distributions, where 9 is the 
angle between a meson and the positron beam. In Figs. 2 (g,h,i) we show that the angular 
distributions of the data events are in excellent agreement with the MC distributions in all 
three cases for ip{2S) — > tt+tt^, K^K~, and KsKl decays. 

The observed peaks for tt+tt^ and K^K~ decays contain contributions from continuum, 
or form factor production e+e^ — )■ 7* — )> 7i^7i^,K^K~ in addition to the resonance con- 
tributions. Continuum contribution in KsKl decays can, however, only arise from higher 
order processes, and is expected to be very small. The continuum contributions have to be 
estimated, and subtracted from the observed peak yields of tt+tt^ and K~^K~ before the 
corresponding branching fractions can be determined. This is particularly challenging for 
the n~^7i~ decays. 
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FIG. 2. Event distributions for ip{2S) -^ tt+tt', K+R-, and KsKl- Panels (a,b,c) show X{h) 
distributions in the extended range X{h). Panels (d,e,f) show enlarged X{h) distributions in the 
region in which we fit the data with MC-determined peak shapes and linear backgrounds. Panels 
(g,h,i), show the corresponding angular distributions. Points with errors represent data; shaded 
histograms represent MC fits to the data. In panels (d,e) the dotted histograms indicate continuum 
contributions. 



The main limitation in the measurement of the interference phase difference angle in 



earlier publications came from the determination of the branching fraction for tt+tt" decay, 
B{ijj{2S) — 7> TT+TT^). As mentioned earlier, the three-gluon strong decay ip{2S) — ?■ tt+tt" is 
forbidden by isospin conservation, and the branching fraction B{iIj{2S) — ?> tt+tt^) is conse- 
quently small. The problem of the intrinsically small branching fraction is compounded by 
the fact that no good-statistics measurements of the continuum tt+tt^ contribution were 
available. As a result all three previous measurements had very large (60 — 100)% errors: 

i3(7r+7r-) X 10^ = 80 ± 50 (DASP [15]), 

= 8.4 ±6.5 (BES [6]), and 
= 8 ±8 (CLEO [8]). 

In the DASP [15] and BES [6] measurements no attempt was made to subtract the continuum 
contribution. In our published paper [8] with 3 million ip{2S) (corresponding to integrated 
luminosity J Cdt = 5.6 pb~ ), only 11 tt+tt^ counts were observed. From these the scaled 
continuum contribution of 7 counts, based on data with e'^e~ integrated luminosity of 20.7 
pb^^ taken oS-%jj{2S) at y/s = 3670 MeV, was subtracted to get A^(7r"''7r~) = 4±4. This led 
to the poor determination of B{iIj{2S) -^ tt+tt-) = (8 ± 8 ± 2) x 10"^ [8]. 

Our present analysis has two big advantages over the old analysis. We now have available 
a much larger data set of 24.5 million ip{2S) corresponding to an integrated e~^e~ luminosity 
of 48 pb~^. Also, we are able to make a much better estimate of the continuum contributions 
in the yields of tt+tt^ and K^K' based on our large-statistics form factor measurements 
with luminosity of 805 pb~^ at v^s = 3772 MeV and 586 pb"^ at v^s = 4170 MeV. 

The widths of V^(3770) and V^(4160) are r(V^(3770)) = (27.3±1.0) MeV and r(V^(4160)) = 
(103±8) MeV, respectively [13]. They are about two orders of magnitude or more larger than 
T{'il){2S)) = (0.304 ± 0.009) MeV, and the estimates of the resonance branching fractions 
for the decays of ■?/'(3770) and ■?/'(4160) to tt+tt", K^K~ , and KgKi range from 1 x 10~^ to 
8 X 10^^. These lead to the conclusion that the resonance contributions in tt+tt^ and K^K~ 
decays of ■?/'(3770) and ■?/'(4160) are less than 0.01% of the total, i.e., the total observed counts 
A^'(tot. obs) are entirely due to continuum or form factor contribution. They can therefore 
be confidently extrapolated to estimate continuum contribution in the measured counts at 
iIj{2S). The extrapolation is done as 

^^ Ar(cont,v/i) C{^)e{^) ^, /^v^V ^2) 

iV'(obs, ^/?) C{^/?)e{V^) \Vs' 



where y/s = M{^{2S)) = 3686 MeV and Vs' = 3772 MeV and 4160 MeV. Here C{y/s, Vs') 
are the luminosities, and e{\^, ys') are efficiencies determined by Monte Carlo (MC) simu- 
lations, and the factor {\fs' / y/s)^ is the conventional extrapolation based on the constancy 
of Q'^F{Q'^) for vector meson form factors [14]. The observed K^K^ and tt+tt^ counts at 
ys' = 3772 MeV and 4170 MeV have statistical and systematic uncertainties [12]. Using 
the counts at 3772 and 4170 MeV in Eq. (2) leads to estimated continuum contributions 
at ^/s = 3686 MeV which are 105.7±3.6±4.7 and 109.2±4.9±4.8 counts respectively for 
kaons, and 41.8±2.2±4.3 and 37.6±3.1±3.9 counts respectively for pions (the ffist errors 
are statistical and the second errors are systematic). We use their averages, taking account 
of the fact that systematic errors are correlated, as 106.9±5.5 counts for kaons, and 40.4±4.6 
counts for pions as our best estimates of continuum contributions at 3686 MeV. 

These contributions are illustrated as dotted histograms in Figs. 2 (d,e). 

In Table I, for ip{2S) decays we list the number A^(obs) of counts observed in the 7T~^7r~, 
K^K^ and KsKl peaks, the number A^(cont) of continuum counts estimated as described 
above, the net signal counts A^(signal) = A^(obs) — A^(cont), the event selection efficiencies 
e, and the branching fractions calculated as 

B{il}{2S) -^ PP) = Ar(signal)/[e x N{^{2S))], (3) 

where the number oiilj{2S) is N{iIj{2S)) = 24.5 x 10^. We note that these branching fractions 
have been obtained without taking account of possible interference between continuum and 
resonance contributions. 

We have considered various sources of systematic uncertainties in our branching fraction 
results. As in Ref. [8], for all three decay channels, ip{2S) — ?■ tt^tt^, K^K^ , and KsKl the 
common uncertainties are ±2% in number of iIj{2S), ±1% per track in track finding, and 
±1% per track in charged particle identification. Uncertainties in trigger efficiency are ±1% 
in TT+TT^ and K^K^, and ±2% in KsK^. The systematic uncertainty in determination of 
the factor C in Eq. (2) comes from the uncertainties in the total integrated luminosity values 
of the data taken at ip{2S) and at ^/s = 3772 MeV and y/s = 4170 MeV, which correspond 
to 1% for each. Thus, we assign 1.4% systematic uncertainty to the value C, determined in 
Eq. (2). 

Variation of the total momentum Epj < 75 MeV requirement by ±15 MeV resulted in 
no statistically significant change in B{ip{2S) — )> tt+tt^), and a ±3.5% change in B{ip{2S) — )■ 
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K^K ), which we assign as systematic uncertainty. For ip{2S) — ?■ tt+tt changing the 
requirement E{CC) > 0.42 GeV by ±10% resulted in ±7% change in i3(^(25) -^ tt+tt-). 
The effect of the implementation of -ft'i-related constraints in KsK^ decay was determined 
as the ratio of fitted peak counts in the spectra in Fig. 2(f)/Fig. 2(c). It was determined 
to be 0.835±0.042. The 5% uncertainty in this determination was assigned as a systematic 
error in B{'4){2S) -^ KsKl). 

The total systematic uncertainties are 8.0%, 5.2%, and 6.5% for ip{2S) — )> -n^-K~ , K~^K~ , 
and KsKl decays, respectively. 

The resulting branching fractions are 

B{i){2S) -^ TT+TT-) = [7.6 ± 2.5(stat) ± 0.6(syst)] x 10-^ 
B{^{2S) -^ K+K-) = [7.48 ± 0.23(stat) ± 0.39(syst)] x 10-^ 
B{ip{2S) -^ KsKl) = [5.28 ± 0.25(stat) ± 0.34(syst)] x 10"^ (4) 

These are listed in Table I. In Table II the errors are listed with the statistical and sys- 
tematic errors combined in quadrature, together with results from previous investigations 
by DASP [15], BES [7, 16], and CLEO [8]. The uncertainties in our branching fractions 
results are factors two to five smaller than those in the published results. In Table II, we 
also list the phase angle difference 6{iIj{2S)) = (113.4 ± 11.5)°, calculated using Eq. (1). All 
the published values of 6{ip{2S)) are also listed, as recalculated using Eq. (1). 

In Sec. V we present the determination of the phase angle differences using a MC method 
which allows us to take account of the distributions in the values of the branching fractions. 

IV. DETERMINATION OF B(J/tp -^ tt+tt", K+K, KsKl) AND 5(J/^) 

The 24.5 million ip{2S) in our data set lead to 8.6 million J/ip events tagged by tt+tt" 
recoil in the decay ip{2S) — )> tt+tt^ J/-*/;. This sample is automatically free of any contami- 
nation of e+e^, fi^fi^, K^K^ , or pp events produced in direct measurement at ^/s = 3097 
MeV. The subsequent J/ip decays to tt+tt^, and K^K~ also do not contain any continuum 
contributions. As such, these data are cleaner and simpler to analyze than the data from 
e^e^ collisions at ^/s = 3097 MeV. 

As stated earlier, to select events for J/ip decays our event selection criteria are modified 
to require four charged particles instead of two. Recoil mass is then constructed for every 
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FIG. 3. For the decay il){2S) — )• -k'^ it~ J / ip , J/if) — )■ KsKl, Ks —5- vr+vr", distributions of the recoil 
invariant mass against vr+vr^ pair which does not come from the Ks decay. Points are for the data, 
shaded histogram corresponds to the signal MC. Normalization is arbitrary. 

pair of two oppositely charged particles. This is dominated by the production of J/ip, as 
shown in Fig. 3 for the channel J/ip — )■ KsK^, which is similar to those obtained for the 
channels J/ip — )■ tt+tt^, K^K~ . We define the clean J /ip sample as consisting of events 
with M(recoil) = M{J/ip) ± 10 MeV. Selection of events for J/ip -^ tt+tt^, K+K', and 
KsKl is done exactly in the same manner as described for iIj{2S). As mentioned earlier, 
the most efficient E{CC) cut to reject ji'^n" from J /ip decays is to require that each pion 
satisfy E[CC) > 350 MeV. The MC-estimated muon contamination with this requirement 
is < 1%. 

Figures 4 (g,h,i), show that the angular distributions of the data events for all thr ee J/ijj 
decays are in excellent agreement with the MC sin^ 9 distributions, as in the case of 'ip{1S) 
decays in Figs. 2 (g,h,i). 

The systematic errors for J /ip decays were determined in exactly the same manner as 
for ipipiS) decays. Their totals are 9%, 6.8%, and 5.5% for tt+tt^, K^K^, and KsKl, 
respectively. 

Figure 4 shows the results for J/ip decays to vr+vr", K~^K~, KsKl as Fig. 2 does for 
the corresponding ip{2S) decays. Figures 4 (a,b,c) show the data in the extended range 
of X{h) with arbitrarily normalized MC predictions. The distribution of the Ks events 
as a function of X{Ks) in the rest frame of J/ijJ, shown in Fig. 4 (c), needs comment. 
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FIG. 4. Event distributions for J/tp — t- vr+vr", K~^K~, and KsK^. Panels (a,b,c) show X{h) dis- 
tributions in the extended range X{h) = 0.8 — 1.1. Panels (d,e,f) show enlarged X{h) distributions 
in the region X{h) = 0.94— 1.06 with fits with MC determined peak shape and linear background. 
Panels (g,h,i) show the corresponding angular distributions. Points with errors represent data, 
shaded histograms represent MC fits to the data. 



The KsKl peak at X{Ks) = 1.0 is clearly separated from the smaller peak at X ~ 0.94 



-*o 



which arises from the J /ip decays into K (892)i^5, K*^ — > Kltt'^, Ks -^ tt+tt , despite 7r° 
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N{ohs) Af(cont) Af(signal) e (%) B{tlj{2S) -^ PP) x 10^ xVdof 



^(25") ^ TT+vr^ 70.8 ±8.8 40.4 ± 4.6 30.4 ± 9.9 16.4 
^-(25") -^ K+K- 1431.3 ± 39.4 106.9 ± 5.5 1324.4 ± 39.8 72.4 
^-(25") -^ KsKl 478.0 ± 23.0 478.0 ± 23.0 53.5 



7.6±2.5±0.6 0.68 

74.8±2.3±3.9 1.11 

52.8±2.5±3.4 1.00 



BiJ/ij) -^ PP) X 10^ 



J/V'^vr+vr- 137.6 ±11.8 

J/V' ^ i^+K" 1057.7 ± 32.8 
J/'ip -^ KsKl 334.3 ± 19.3 



137.6 ±11.8 10.9 
1057.7 ±32.8 43.1 
334.3 ±19.3 21.5 



1.47±0.13±0.13 1.09 

2.86±0.09±0.19 1.00 

2.62±0.15±0.14 1.03 



TABLE I. Fit results for ip{2S), J/ijj — )■ tt+tt , K^K , and KsKl decays, and the corresponding 
branching fractions. 

rejection. The clear separation of the KsK^ peak from the K (892)i^5 peak is confirmed by 
the MC simulation whose result is superposed on the data in Fig. 4 (c). In Figs. 4 (d,e,f), 
the event distributions of Figs. (a,b,c) are shown in detail in the smaller region of X{h). 
Also shown are the fits using MC-determined peak shapes and linear backgrounds. The 
fits give 137.6 ± 11.8, 1057.0 ± 32.8, and 334.3 ± 19.3 counts for tt+tt", K+K^, and KsKl, 
respectively. These compare with 84, 107 and 74 counts in the Mark III measurements with 
a factor three smaller sample of J/ip [10]. These counts, the MC determined efficiencies e, 
and N{J/%1)) = (8.57 ± 0.07) x 10^ [17], lead to the branching fractions listed in Table I. 
Thus the final branching fractions are 

B{J/^ -^ TT+TT-) = [1.47 ± 0.13(stat) ± 0.13(syst)] x 10-^ 
B{J/^p -^ K+K-) = [2.86 ± 0.09(stat) ± 0.19(syst)] x 10-^ 
B{J/^ -> KsKl) = [2.62 ± 0.15(stat) ± 0.14(syst)] x 10"! (5) 

These are listed in Table I, and in Table II with the statistical and systematic errors combined 
in quadrature. The above branching fractions for J/ip — )> tt+tt^ and K^K~ are consistent 
with those reported by Mark III [10], and have factors two and three smaller uncertainties, 
respectively. Our branching fraction for J /ip — )■ KsKl based on 334±19 well resolved events, 
as shown in Fig. 4 (f), is a factor 2.6 larger than B{J /i) -^ KsKl) = (1.01 ± 0.18) x 10"^ 
reported by Mark III with 74 identified counts, obtained by "stringent cuts" to remove 
J/ip -^ p°7r° and J/ip -> KsK* {898) decays, and is ~ 40% larger than B{J/iIj -^ KsKl) = 



12 



DASP [15] BES [6, 7] CLEO [ 
1979 2004 2005 



This analysis This analysis 
MC result 



^(V'(2S') ^ vr+vr-) X 10^ 80 ± 50 8.4 ± 6.5 8±8 
l3{'tp{2S) ^K+K-)xW^ 100 ±70 61 ± 21 63 ± 7 
Bii;{2S) -^ KsKl) X 10^ — 52.4 ± 6.7 58.0 ± 9.0 



,5(^(25)) 



7.6 ±2.6 

74.8 ±4.5 
52.8 ±4.2 



(91 ±35)° (87 ±20)° (113.4 ±11.5)° (110.5+g^5°)° 



Mark III[10] BES [16] 



This analysis 



This analysis 
MC result 



B{J/%1) -^ 7r+^-) X 10^ 
^(J/V' -^ K+K-) X 10^ 
B{J/i; -^ KsKl) X 10^ 



1.58 ±0.25 
2.39 ±0.33 
1.01 ±0.18 

(88 ±11)° 



1.82 ±0.13 



1.47 ±0.18 

2.86 ±0.21 

2.62 ±0.21 

(73.6 ±5.6)° 



+5.0x0 



(73.5^^:^) 



TABLE II. Summary of results for ip{2S) and J/ip decays to pseudoscalar pairs: branching fractions 
and the phase difference angles 6{iIj{2S)) and 6{J/iIj) using central values of branching fractions. 
The BES [7] and CLEO [8] results for 6{tp(2S)) have been recalculated to correspond to the internal 
angle of the amplitude triangle. In the last column S^ip) results based on Monte Carlo calculation 
described in the text are presented. 



(1.82 ± 0.13) X 10"^ reported by BES II [7] with 2155 ± 45 identified events. In the large 
statistics BES II measurements the peak due to J/^p — )■ K {892)Ks, K — )■ Kltt^, Ks — >■ 
7r^7r+, overlapped with the direct J /ip — )■ KsKl, Ks — )> tt+tt^ peak, and strong cuts had 
to be made to resolve the two peaks. As shown in Figs. 4 (c,f), in our measurements 
the two peaks are completely resolved. Further, MC calculations confirm that with our 
selections the decay J /ip — )■ vr^p*^, p^ — > tt+tt" also does not make any contribution under 
the J/i) -^ KsKl, Kg -^ Tr+Tr" peak at X{Ks) = 1.0. 

As for iIj{2S), we calculate S{J/ip) using Eq. (1), and obtain S{J/ip) = (73.6 ± 5.6)°, as 
compared to (88 ± 11)° obtained using the branching fractions measured by Mark III. These 
are listed in Table II. 
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FIG. 5. (Left) The relative phase angle distributions Si^ip) for ip{2S) and J/ip obtained from toy MC 



simulations. The results are (5(V^(25)) = (llO.Sl^^g'^)" and SiJ/ip) = (73.5 



:;+5.0^o 
-4.5/ • 



(Right) Toy MC 



distribution of the difference of relative phase angles for ip{2S) and J/ip, A5 = 6{ip{2S)) — 6{J /ip) = 

(sr.otllir. 

V. MONTE CARLO BASED EVALUATION OF 6{ij) AND THE DIFFERENCE 

A(5 = 6{'ilj{2S)) - 5{J/^l>) 

The individual measured branching fraction values have distributions which are conven- 
tionally stated in terms of la errors, as listed in Table I. Using only the central values to 
evaluate 5{'4)) according to Eq. (1) is not correct. The more correct procedure is to make 
Monte Carlo evaluation of Eq. (1) taking account of random associations of the branching 
fraction values in the distributions for the three decays. We have made such toy MC eval- 
uations of both 5{il){2S)) and 6{J/iIj). As shown in Fig. 5 (left), the large error (±30%) in 
ip{2S) — 7> TT+TT^ branching fraction leads to a very asymmetric MC distribution for 6{iIj{2S)), 
whereas the much smaller error (±12%) for J/ip — )> tt+tt" results in a much smaller asym- 
metry in the distribution for 6{J/iIj). If we adopt the usual definition of the la error as that 
which includes 68% of the area on each side of the peak of a distribution, our results are 



5{i;{2S))pp = (110.5™)°, and 5{J/iP)pp = (73.5+^°)^ 



(6) 
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The difference, whose MC distribution is illustrated in Fig. 5 (right), is 

A6 = 6{ij{2S))pp - 6{J/ij)pp = {37.0t\liy. (7) 

We consider the above estimates of 6{ip{2S))pp, 6{J/ip)pp, and their difference to be our 
final results. 

VI. CONCLUSIONS 

We have made large statistics measurements of the branching fractions for the decays of 
ip{2S) and J/ip to pseudoscalar pairs tt+tt^, K^K^, and KsKl. Our branching fraction 
results have errors which are factors two to five smaller than the previously published results. 
Using these branching fractions we have made calculations of the phase angle differences 
between the electromagnetic and strong decays for both ip{2S) and J/ip, taking proper 
account of the distributions of the branching fraction values. Our results are nearly a factor 
two more precise than the previously published results. 
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